Determination of the Equation of State of a Polarized Fermi Gas at Unitarity 
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We report on the measurement of the equation of state of a two-component Fermi gas of ^Li 
atoms with resonant interactions. By analyzing the in situ density distributions of a population- 
imbalanced Fermi mixture reported in the recent experiment [Y. Shin et al, Nature 451, 689 (2008)], 
we determine the energy density of a resonantly interacting Fermi gas as a function of the densities 
of the two components. We present a method to determine the equation of state directly from the 
shape of the trapped cloud, where the fully-polarized, non-interacting ideal Fermi gas in the outer 
region provides the absolute calibration of particle density. From the density profiles obtained at 
the lowest temperature, we estimate the zero-temperature equation of state. 
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Interacting fermions are a paradigm of modern physics, 
encompassing superconductivity and superfluidity. One 
interesting situation arises when the constituents interact 
resonantly, i.e., the scattering length for the free fermions 
diverges. At this so-called unitarity limit, the behavior of 
the system becomes universal, being independent of the 
nature of the interactions. Ultracold atomic Fermi gases 
near Feshbach coUisional resonance present a model sys- 
tem for studying strongly interacting fermions [1]. Re- 
cently, the phase diagram of a two-component Fermi gas 
with resonant interactions has been experimentally es- 
tablished [5], showing that at zero temperature the sys- 
tem undergoes a first-order quantum phase transition 
from a fully-paired superfluid to a partially-polarized nor- 
mal gas when the imbalance between the two spin compo- 
nents exceeds a critical value, called the Chandrasekhar- 
Clogston limit of superfluidity [31 H] . 

The nature of the partially-polarized normal phase, 
however, is still a subject of investigation. The spectral 
shift observed in the minority rf excitation spectrum has 
been interpreted as the existence of 'pairing' in the nor- 
mal phase [S] , but several theoretical studies for a highly 
polarized system, e.g. a single minority atom in a major- 
ity Fermi sea, suggest that the system is well-described 
as a normal Fermi liquid, where the minority atoms are 
associated with weakly interacting quasiparticles [SI El [5] • 
This picture seems to be supported by the experimental 
observation that the shape of the minority cloud in the 
normal phase is similar to a free Fermi gas [21 [9j . It has 
been speculated that exotic pairing states might exist in 
the partially-polarized phase [TP] . 

In this Letter, we determine the equation of state of 
a polarized Fermi gas at unitarity from the the in situ 
density proflles of a population-imbalanccd Fermi mix- 
ture confined in a harmonic trap. Since the variation of 
the external trapping potential across the sample scans 
the chemical potential, in principle, the density informa- 
tion of a single sample contains the whole equation of 
state [m 112]. We present a method to determine the 



equation of state directly from the shape of the trapped 
cloud. Because of its exactly-known thermal properties, 
a fully-polarized, non-interacting ideal Fermi gas in the 
outer region provides the absolute density calibration. 
The equation of state of a polarized Fermi gas can be 
parameterized using a normal Fermi liquid description, 
which includes the binding energy of a single minority 
atom resonantly interacting with a majority Fermi sea, 
the effective mass of the quasiparticles, and its correc- 
tion. This work is the first quantitative study of the 
thermodynamic properties of the polarized normal state 
with strongly interactions, finding reasonable agreement 
with recent calculations [^ ITIIF HTTI [H] . 

For infinite scattering length, the unitarity limit im- 
plies that all interaction energies scale with the Fermi 
energies of the components epi — /2m{6TT^ni)'^/^ |13j . 
where h is the Planck constant divided by 27r, m is the 
particle mass, Ui is the density of component i, and 
i = 1,2. As a result, a simple dimensional scaling ar- 
gument implies that the energy density f(ni,n2) of a 
two-component Fermi gas can be parameterized as 

£{m,n2)^^a[mg{x)]''/\ (1) 

introducing a dimensionless universal function g{x) |12j . 
where a = (67r^)^/^?i^/2TO and x = n-ilrix is the density 
ratio of the two components. Without loss of generality, 
< x < 1 due to the symmetry of the two components. 
From the chemical potential relation [ii ~ d£ /drii, the 
universal function can be expressed as 

where y — fJ-2/li-i is the chemical potential ratio of the 
two components. 

The main result of this paper is the measurement of 
the universal function g(x) for a resonantly interacting 
Fermi gas. When a Fermi mixture is confined in a har- 
monic trap, V{r) oc r^, the local chemical potential is 
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given as Hi{r) = /Xio — V(r), where fiio is the global chem- 
ical potential with respect to the bottom of the trapping 
potential. The global chemical potential of the majority 
component /xig is determined from the radius of the ma- 
jority cloud i.e., fiio = V{Ri). Then, the chemical 
potential ratio y{r) — iJL2{r) / ^,i{r) is given as 
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where yo — M2o/a'io- With population imbalance, i.e. 
yo < 1, x{r) and y{r) vary over the sample. From Eq. (|2]), 
the spatial correlation of the local particle densities and 
the local chemical potentials in a trapped sample deter- 
mines g{x) under the local density approximation. 

We estimate the zero-temperature equation of state 
by analyzing the density profiles obtained at the lowest 
temperature (T/Tp « 0.03 where Tp is the Fermi tem- 
perature of the majority component) in Ref. 2 . A spin 
mixture of the two lowest hyperfine states of ^Li atoms 
was prepared in a 3D harmonic trap on a broad Fesh- 
bach resonance, located at a magnetic field of 834 G [T3], 
resonantly enhancing the interactions between the two 
spin states. The detailed description of the experimen- 
tal procedure for the sample preparation and the signal 
processing has been provided in Ref. [2 . 

Figure [T] shows the spatial structure of a resonantly in- 
teracting Fermi mixture in a harmonic trap. According 
to the zero-temperature phase diagram [2\ , three distinc- 
tive spatial regions can be identified in the inhomoge- 
neous sample with population imbalance (j/g < 1). When 
the chemical potential ratio at center j/q is larger than a 
critical value j/c for the superfluid-to-normal phase tran- 
sition, the sample has (I) a fully-paired superfluid core 
{x = 1, y > yc) surrounded by (II) a partially-polarized 
normal gas (0 < a; < Xc, y < yc), showing a discontinuity 
in the density ratio x at the phase boundary. In the outer 
region where the minority component is completely de- 
pleted, (III) a fully-polarized non-interacting Fermi gas 
forms (x — 0, y < y^ < yc)- The radii, Rc and i?2 are 
the I-II and II-III transition points, defining y{Rc) = yc 
and y(i?2) = Vm, respectively. 

The non-interacting ideal Fermi gas in the outer (r > 
R2) region provides a reliable method to measure the lo- 
cal chemical potential /ii(r) in the strongly interacting, 
inner (r < R2) region. Since epi = /^i in the outer region, 
the extension of the non-interacting ideal Fermi distribu- 
tion no(r), fit to the outer region, into the inner region 
gives the local chemical potential as fJ-i{r) — an^^^, con- 
sequently ^i/epi = (no/ni)^/^. This method allows to 
measure the equation of state directly from the shape of 
the cloud without any absolute calibration for particle 
density. Furthermore, when a sample has a superfluid 
core, i.e. the whole range of the density ratio < x < 1, 
a single shot image of the sample can provide all infor- 
mation for the determination of the equation of state. 
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FIG. 1: (Color online) Spatial profiles of a population- 
imbalanced Fermi mixture confined in a harmonic trap at 
unitarity. (a) The averaged column density profiles and (b) 
the reconstructed 3D profiles at the lowest temperature from 
Ref. ;2: (green (gray): majority, blue (dark gray); minority, 
black: difference). Tii and R2 are the radii of the majority 
(spin 1) and the minority (spin 2) cloud, respectively. The 
core radius Rc was determined as the kink and peak position 
in the column density difference. The black dotted line in (a) 
is a zero-temperature Thomas-Fermi (TP) distribution fit to 
the outer region (r > R2) of the majority column density pro- 
file and the black dotted line in (b) is the corresponding 3D 
distribution, no- The sample has three regions: (I) The su- 
perfluid core region (0 < r < Rc), (II) the partially-polarized, 
intermediate region {Rc < r < R2), and (III) a fully-polarized 
non-interacting outer region {R2 < r < Ri). (c) The density 
ratio of the two components, x{r) = n2/ni. (d) The majority 
density normalized by the reference density, tiq. 



In the experiments, the total population imbalance was 
controlled to be 5 = (iVi - A^2)/(iVi + N2) = 44(4)% 
less than the critical imbalance i5c(« 75%) [15j \W\ to 
have a superfluid core, i.e. yg > yc, where Ni is the 
total atom number of component i. The phase boundary 
Rc, located by the kink and peak position in the column 
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FIG. 2; (Color online) Thermodynamic potential at unitar- 
ity. The universal function g{x) for energy density was con- 
structed as a function of the density ratio x of the two com- 
ponents, from five independent measurements of the density 
profiles (Fig. The critical chemical potential ratio and 
the critical density ratio were measured to be yc = 0.03(2) 
and Xc — 0.53(5), respectively (see text for determination 
method). The red (dark gray) solid line is obtained by 
fitting the model equation Eq. (j4| with {^,yc,ym, Xc} = 
{0.42, 0.03, -0.58, 0.53} to the normal region {x < Xc) of g{x). 
A red (dark gray) dashed line in the phase separation region 
{xc < a; < 1) connects the two points {x = 1, g{l) = (2^)^''^) 
and {x = Xc, gixc))- 



density difference profile, was measured to be Rc/Ri = 
0.430(3) and the critical density ratio was measured to 
be Xc = 0.53(5). The reference density no{r) and the 
radius i?i(i?2) were determined from the fit of the outer 
region, r > R2{r > Re) of the majority (minority) column 
density profile to a zero-temperature Thomas-Fermi (TF) 
distribution. 

The well-known zero-temperature thermodynamics of 
a balanced superfluid and a non-interacting Fermi gas 
provide physical constraints on the form of the universal 
function. Since the chemical potential of a fully-paired 
{x — 1) superfluid is proportional to the Fermi energy, i.e. 
IJLs = (mi + M2)/2 = Cefi, we have Hi/epi = 2^/(1 -I- y) 
and g{l) = (2^)'^/^. On the other hand, a fully-polarized 
{x = 0) non-interacting Fermi gas has /xi = epi so 
that g{Q) — 1. The universal parameter ^ can be inde- 
pendently determined from the majority profile by com- 
paring the curvature of the Fermi energy distribution 
£Fi{r) oc ni{r)'^/^ in the core region and in the outer 

region, as ^ = {d?£Fi/ dr'^)r>R^/ {dPepi/ dr'^)r<Ra- This 
determination is, however, technically limited due to the 
low signal-to-noise ratio. In the following analysis we use 
the theoretically predicted value £,th = 0.42(1) [T71 118j. 
confirmed in previous measurements [13 [201 HH Ull 

EH 123 Eg. 

Figure |2] displays the universal function g{x) con- 
structed from the density profiles at the lowest tem- 
perature. The critical value yc was determined to be 
yc = 0.03(2) such that the average value of g{x) for 
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FIG. 3: (Color online) Comparison between the experimen- 
tal data and the model, (a) The 3D density profiles and (b) 
the column density profiles for 5 — 44% (dot-dashed: major- 
ity, dashed: minority, dotted: difference) generated from the 
model equation of state (the red (dark gray) line in Fig. [2]| . 
Experimental data in Fig. [l] are plotted together for compar- 
ison (same color usage). 



X > 0.9 gives {2^th)^/^ = 0.90. This critical value yc 
has been discussed in Ref. 2 to demonstrate the stabil- 
ity of a fully-paired superfluid state at zero temperature. 
At zero temperature, g{x) is not deflned for a homoge- 
neous system with Xc < x < 1. The sparse population 
of the data points in the region of < a; < 1 indicates 
the phase separation in the sample, associated with the 
first-order phase transition. 

The other critical value ym represents the binding en- 
ergy Eb of a single minority resonantly interacting with 
a majority Fermi sea as Eh = lim2,^o+ M2 = ym^Fi- 
y = 9'{x)/{g{x) - xg'{x)) and 5(0) = 1 give y„ = 
^'(0). By fitting 1 + y„iX to g{x) for x < 0.1, we esti- 
mated ym = —0.58(5), which is in good agreement with 
the recent theoretical results of —0.6 [11], —0.58(1) [6], 
-0.6066 [7], -0.54(4) [12], and -0.618 [8]. From the 
definitions, yc = y{Rc) and = 2/(^2)1 the estimated 
critical values {j/cj/m} = {0.03(2), —0.58(5)} suggest 
R2 = 0.707(20)i?i for Rc = 0.43i?i. However, the mi- 
nority radius was measured to be R2/R1 = 0.73(1), sug- 
gesting ym = —0.69(8). We attribute this discrepancy to 
the fact that the shape of the minority density profile in 
the intermediate region cannot be completely captured 
with a zero-temperature TF distribution [27i . 

Following a normal Fermi liquid description [Sj [Tj H] , 
we consider a model for the equation of state of a partially 
polarized Fermi gas in the following form, 

g{xf'^ = 1 + r^VraX + '^^^^x'' for x < Xc- (4) 
3 m* 

satisfying the boundary conditions g{0) ~ 1 and g'{0) = 
ym at a; = 0. The second term corresponds to the 
momentum-independent binding energy for the minor- 
ity atoms and the third term describes the deviation 
from the free particle behavior, regarding m* as the effec- 
tive mass and c as its correction (for the non-interacting 
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case, go{x)^^^ = 1 + x^^^). The equilibrium condition 
for the coexistence of two spatially separate phases re- 
quires that the two phases have the same chemical po- 
tential and pressure at the critical point |28j . imposing 
the boundary conditions at x — Xc such as g{xc) — 
[(1 + Xcyc)/il + VcMl) and g'{x,) = [y,/{l + ycMl)- 
Then, for given values {^,yc,ym,Xc}, this model relies 
on only one free parameter. With {£,,yc,ym, Xc} — 
{0.42,0.03,-0.58,0.53}, the fit of Eq. (|t| to the inter- 
mediate region (x < Xf.) gives 7 = 1.60(13), having 
m* = 1.06 and c = —0.019. Quantum Monte-Carlo 
calculations for small x predicts m* = 1.04(3) with 
7 = 5/3 [IT, which is very close to the observed behavior. 
Figure [3] displays the density profiles generated from the 
model equation of state, together with the experimental 
data. 

Our observation of 7 « 5/3 and the small change in the 
effective mass suggests that a polarized Fermi gas with 
resonant interactions can be described as a normal Fermi 
liquid with weakly interacting quasiparticles. However, it 
is an open question whether the Fermi liquid description 
is still valid for high minority concentrations, where the 
Pauli blocking effect of the minority Fermi sea might play 
an important role. The possibility of the exotic ground 
state of a partially polarized system has been suggested 
by the recent observation of the temperature-dependent 
spectral shift in the minority rf excitation spectrum [5] . 
We note that it is not clear how to distinguish possible 
exotic states fW via the equation of state. More experi- 
mental studies for microscopic properties of the system, 
e.g. majority rf spectroscopy, are necessary to clarify the 
issue. 

In conclusion, we measure the equation of state of a 
two-component Fermi gas with resonant interactions by 
analyzing the in situ density distributions of the trapped 
sample. In a similar way, the density profiles at finite 
temperature may reveal the excitation spectrum of the 
system [53] • 
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